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ABSTRACT 

Aims. It is commonly assumed that the massive winds of AGB stars are dust-driven and pulsation-enhanced. However, detailed 
frequency-dependent dynamical models which can explain the observed magnitudes of mass loss rates and outflow velocities have 
been published so far only for C-stars. This letter reports on first results of similar models for oxygen-rich AGB stars. The aim is to 
provide a better understanding of the wind driving mechanism, the dust condensation sequence, and the role of pulsations. 
Methods. New dynamical models for dust-driven winds of oxygen-rich AGB stars are presented which include frequency-dependent 
Monte Carlo radiative transfer by means of a sparse opacity distribution technique and a time-dependent treatment of the nucleation, 
growth and evaporation of inhomogeneous dust grains composed of a mixture of Mg2Si04, SiC>2, ALO3, TiC"2, and solid Fe. 
Results. The frequency-dependent treatment of radiative transfer reveals that the gas is cold close to the star (700 - 900 K at 1.5-2/?*) 
which facilitates the nucleation process. The dust temperatures are strongly material-dependent, with differences as large as 1000 K for 
different pure materials, which has an important influence on the dust formation sequence. Two dust layers are formed in the dynamical 
models: almost pure glassy AI2O3 close to the star (ri 1.5/?*) and the more opaque Fe-poor Mg-Fe-silicates further out. Solid Fe or 
Fe-rich silicates are found to be the only condensates that can efficiently absorb the stellar light in the near IR. Consequently, they play 
a crucial role for the wind driving mechanism and act as thermostat. Only small amounts of Fe can be incorporated into the grains, 
because otherwise the grains get too hot. Thus, the models reveal almost no mass loss, and no dust shells. 

Conclusions. The observed dust sequence ALO3 — » Fe-poor Mg-Fe-silicates for oxygen-rich AGB stars having low — » high mass 
loss rates is in agreement with the presented model and can be understood as follows: ALO3 is present in the extended atmosphere of 
the star below the wind acceleration region, also without mass loss. The Mg-Fe-silicates form further out and, therefore, their amount 
depends on the mass loss rate. The driving mechanism of oxygen-rich AGB stars is still an unsolved puzzle. 

Key words. Hydrodynamics - Radiative transfer - Stars: winds, outflows - Stars: mass loss - Stars: AGB and post-AGB 



■ 1. Introduction In contrast, the a posteori frequency-dependent radiative 

transfer analysis of non-linear pulsation models with simplified 

, , . . . , , . . dust formation theory by Ireland & Scholz (2006) did not find 

The mass loss mechanism of AGB stars and red supergiants is , .. , , ,,, n \ . . „ C ( ~. , 

, ,. ■ • , , , T 1 11 much radiation pressure on dust (ALO3 and Mg2j-Fe2-2j:Si04) 

a long-standing astrophysical problem. In the carbon-rich case, „ . • , , ,• .- , ,, 

,. , 3 . , , in O-nch Mira variables, with radiative accelerations as small as 
an extraordinary condensate exists (amorphous carbon) which 



is very stable, i. e. it can exist already close the star, and is 
very opaque in the optical and near IR spectral region. Detailed 
dynamical models with time-dependent dust formation theory 



0.08 to 0.29 times the gravitational deceleration. 

Recent mid-IR observations of O-rich AGB stars in globu- 
lar clusters with Spitzer (Lebzelter et al. 2006) and of galactic 



, I,,,--.- , bulge AGB stars with ISO (Blommaert et al. 2006) show a clear 

(Winters et al. 2000) show that the formation of amorphous , . , , , . , „ , , , 

, -j rc-^j - j- correlation between the kind of condensate and the mass loss 

carbon can provide sufficient radiation pressure to drive mas- „ , „ , , , , 

„ r . , . T . . „„ . rate M, called the observational dust condensation sequence • 

sive outflows, consistent with the basic characteristics of C-star . , , ,-. . . , . , „ , .1 , . , 



winds. This result has been confirmed by dynamical models with 
frequency-dependent radiative transfer by Hofner et al. (2003). 



stars with low M show mainly AI2O3, whereas stars with higher 
M show increasing amounts of Mg-Fe-silicates. From Midi in- 
terferometry of the red supergiant a Ori, Verhoelst et al. (2006) 



j j . , 11. concluded that ALO3 grains are already present at radial dis- 

However, in the oxygen-rich case, no such condensate ex- „ , % _ 3 r 

• * tu * * ui Ti a rt A. n t T-i tances as small as r- 1.5k*. 

ists. The most stable metal-oxides like AI2O3 are too rare. The 

abundant pure silicates like Mg2Si04 are already less stable and 

almost completely transparent around 1 fim where most of the 2 Jhe model 

stellar flux escapes. Solid Fe and Mg-Fe-silicates are opaque but 

even less stable. Stationary models of dust-driven O-rich AGB Hydrodynamics is solved by using the FLASH-solver (Fryxell 

star winds with grey radiative transfer (Ferrarotti & Gail 2006) et al. 2000) in spherical symmetry, including gravity and self- 

and dynamical models with pulsation and grey radiative transfer developed modules for radiation pressure on dust & molecules 

(Jeong et al. 2003) nevertheless came to the conclusion that the and radiative heating /cooling (see Woitke 2006a for details), 

winds of O-rich AGB stars are dust-driven, where the stellar pul- We use here, however, a piston approximation as inner boundary 

sation helps to provide the necessary density conditions to form condition to simulate the pulsation of the star, and a new equa- 

the dust close to the star ("pulsation-enhanced"). tion of state for a mixture of H + , e~, H, H2, He and other atomic 
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Fig. 1. Extinction efficiencies over particle radius Q ext /a in the 
Rayleigh-limit of Mie theory according to the Jena optical data 
base. The data is partly log-log extrapolated. 

metals in LTE, including ionisation and dissociation potentials, 
and vibrational and rotational excitation energies of H2. 

Radiative transfer: For the models presented in this letter, we 
have developed a new frequency-dependent Monte Carlo radia- 
tive transfer technique which allows for arbitrarily high opti- 
cal depths (Woitke 2006b). In frequency-dependent stellar at- 
mospheres, the gas is always optically thick at least in some 
wavelengths (~ 10 5 ), which is a big problem for standard MC 
techniques. The method is coupled to the hydrodynamics and 
can be used also for 2D models. 

Opacities: The basis for our radiative transfer treatment 
are monochromatic molecular gas opacities from the Marcs 
stellar atmosphere code (J0rgensen et al. 1992), extracted by 
Helling et al. (2000). The frequency space is subdivided into 
five spectral bands with two opacity distribution points in each 
band, resulting in altogether 5x2 effective wavelengths sampling 
points (for the traditional ODF approximation see e.g. Carbon 
1979). High and low mean opacity values are tabulated for each 
spectral band during the initialisation phase of the program in 
such a way that they simultaneously result in the correct Planck 
and Rosseland band-mean gas opacities. The details will be ex- 
plained in another paper (Woitke 2006b). 

Dust opacities are calculated in the Rayleigh limit of Mie 
theory according to the Jena optical data base, kindly provided 
by Th. Posch (see Fig. [Q. The total dust extinction coefficient 
[cm 2 /g] is assumed to be given by 1 



Fig. 2. Hydrostatic initial model (full) in comparison to a spher- 
ical Marcs model (dashed) and a grey model (grey r g -line) for 
M+ = 1 M Q , = 2800 K, logg = -0.6 (L* = 6048 L ), Z = 1. 
^Ti02 shows the supersaturation ratio of Ti02, indicating that nu- 
cleation is already possible very close to the star. Long-dashed 
graphs belong to the r.h.s. axis. 

oxygen-rich condensates have a "glassy" character. They are al- 
most transparent in the optical and near IR but opaque in the mid 
IR where the strong vibrational resonances are situated. In con- 
trast, solid Fe and Fe-rich silicates are opaque even in the optical. 
The monochromatic dust opacities are subject to the same aver- 
aging procedure to result in high/low band-mean dust opacities 
as described above for the gas opacities. 

Dust formation is described by a system of differential mo- 
ment equation explained in (Helling & Woitke 2006) consider- 
ing the growth and evaporation of inhomogeneous dust grains 
composed of a mixture of Mg2SiOi, SiC<2, AI2O3, TiC»2, and 
solid Fe (13 growth/evaporation reactions). The nucleation rate 
of (TiC<2)N clusters is adopted from Jeong (2000). The molecu- 
lar concentrations entering into the calculation of the nucleation 
and growth rates are calculated by a small neutral equilibrium 
chemistry for 1 1 atoms (H, He, C, O, N, Mg, Al, Si, S, Ti, Fe) 
and 33 molecules (H2, CO, C0 2 , OH, H z O, CH 4 , N 2 , CN, HCN, 
NH,, H 2 S, SiS, SO, HS, SiO, SiH, SiH 4 , Si0 2 , SiN, S0 2 , MgH, 
MgS, MgO, MgOH, Mg(OH) 2 , FeO, Fe(OH) 2 , AlOH, A10 2 H, 
A1 2 0, A1H, TiO, TiQ 2 ). 



-dust 
K A, ext 



(l) 3. The static solution 



where L3 is the third moment of the dust size distribution func- 
tion and Vj/Vtot 1S the volume fraction of solid material s in 
the dust component. The extinction efficiencies over particle ra- 
dius of the pure solids s are shown in Fig. \l\ Note that most 

1 According to our assumption of dirty grains, an application of ef- 
fective medium theory would be more appropriate, which will be exam- 
ined in a future paper. First results show that the effective extinction is 
stronger than the simple volume-means used in this paper. 



The hydrostatic, dust-free solution (see Woitke 2006a for equa- 
tions and further details) used as initial model for the dy- 
namical models is shown in Fig. |2] The obtained degree of 
agreement with the Marcs model is quite remarkable for such 
a rough 10-wavelengths-point treatment of radiative transfer. 
Similar results have been obtained by Hofner et al. (2003) and 
Helling & J0rgensen (1998), using a sparse opacity sampling 
technique with about 100 sampling points. The gas tempera- 
ture structure T g (r) shows a pronounced "step" of almost 1000K 
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around 1.25 7?* (the Marcs model is actually not extended 
enough to reveal this step completely). It is clear, however, that 
this is not just a "surface effect". Outside of this r g -step even the 
strongest molecular lines become optically thin (in the hydro- 
static case) and the line blanketing effect works at full strength. 
The grey model fails completely in predicting this step which is 
very meaningful for the dust formation. The agreement with the 
Marcs results concerning the mass density p(r), the gas pressure 
p(r), the mean molecular weight fi(r), and the acceleration by ra- 
diation pressure on molecules divided by gravity r gas = a m &lg(r) 
is also satisfying. 

4. Rough estimates of the dust acceleration 

The calculation of the spectral mean intensities J \ (r) and spec- 
tral fluxes (r) in the static model (see Fig. [5) allows for a quick 
estimate of the maximum possible radiative acceleration by dust 
if the dust is still optically thin. Table^shows the resulting dust 
temperatures T d of several pure condensates 



Table 1. Calculated dust temperatures T d (first row) and dust 



■■?ZBAT d )dA 



and the radiative acceleration by dust divided by gravity 



rdustM 



GM(r) 



(2) 



(3) 



at three selected distances from the star. For each condensate 
we take the maximum possible dust volume per mass £3, max 
given by element conservation constraints, e.g. Z^^° 4 — 
Min{±e Mg , e si , |e } V Mg2 sio 4 "<h>/p where «<h> is the hydrogen 
nuclei density, the abundance of element k and VMg 2 sio 4 the 
monomer volume of Mg2SiC<4 (see Helling & Woitke 2006). 

The results shown in Tabled demonstrate that the dust tem- 
peratures 7d are strongly material-dependent, with differences as 
large as 1000 K at the same distance from the star, which is a re- 
markable result. All condensates (except solid Fe) have strongly 
peaked mid-IR resonances which are situated just around the 
maximum of the local Planck function - they work perfect for 
radiative cooling. In contrast, the glassy character of the oxides 
and pure silicates like AI2O3, S1O2, Mg 2 SiC>4 and MgSiC"3 (the 
low absorption efficiencies at optical and near-IR wavelengths, 
see Fig. ^ prevent efficient heating by the star. Consequently, 
the pure glassy condensates can exist astonishingly close to the 
star (see also Woitke 1999). 

For the same reasons, radiative pressure on all glassy con- 
densates is negligible! It is without effect for the wind acceler- 
ation mechanism whether for example Mg2SiC»4 condenses out 
or not. The only dust species that can potentially drive a stellar 
outflow are solid Fe and Fe-rich silicates like MgFeSiC<4. 

The unavoidable consequence of the spectral characteristics 
of oxygen-rich dust is that radiative acceleration must be paid 
for by radiative heating, i. e. dust species capable of driving a 
stellar wind (TS 1) are hot and can only exist at a relatively large 
distance from the star (e. g. r ^5 marked in bold in Table[Q. 
For comparison, amorphous carbon is so opaque and stable (in 
a C-rich gas) that it could accelerate the gas already from 2 
onwards, with 20 x the local gravity in this model. 

5. Results of the dynamical models 

The first results of the dynamical models showed almost no mass 
loss (M ^ 10~ 10 M /yr), just some erratic large-scale and long- 
term excursions for which the mass loss rate is actually difficult 



radiative accelerations Tdust = a< zf/g (second row) in case of 
full condensation of several condensates into small particles in 
the static model (see Fig. |2j. The resulting dust-to-gas ratio 
Pdust/Pgas is shown in the middle column. Temperatures values 
with * mark thermally unstable condensates. 



solid material 
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to measure. We then approached rather extreme stellar param- 
eters (Mi, = 1 M Q , T+ = 2500 K, L+ = 10000 L@), still without 
success. Finally, in order to see how a dust-driven wind could 
look like, we arbitrarily enhanced the radiative acceleration by 



Fgas "t" 5 Fdust 



(4) 



The results of this simulation is shown in Fig. [3] 

The models (also those without enhanced F) show extended 
warm molecular layers, truncated by the r g -step described in 
Sect.[3]due to the line blanketing effect. The pulsation of the star 
leads to a time-dependent extension of these layers to roughly 
1.5 — 2R*. If dust forms, it fills in the gas opacity gaps in fre- 
quency space which reduces the line blanketing effect. 

The formation of seed particles (see the nucleation rate 
J+/ft(H) in Fig. [3} happens right above these molecular layers, 
i. e. very close to the star. The gas is cold here [T g « 700 -900 K) 
which is not revealed by models using grey radiative transfer 
(Jeong et al. 2003, Ferrarotti & Gail 2006). 

According to the model, two dust layers develop: almost pure 
glassy AI2O3 grains close to the star (r ^ 1 .5 even partly in- 
side the warm molecular layers up to T g « T d w 1500 K), and 
further out the more opaque Mg-Fe-silicates grains which grow 
on top of the AI2O3 particles with a very small, steadily increas- 
ing Fe content. 

The temperature of the dirty dust grains is controlled by 
the iron content, which has already been noted by Tielens et 
al. (1998). The volume fraction of solid iron inclusions relaxes 
quickly to a level where a further increase would cause too much 
radiative heating and thus thermal re-evaporation of the solid 
iron inclusions 



(5) 



where r| e is the sublimation temperature of solid iron. This is 
a robust and recurrent result concerning many simulations, be- 
cause there is a stable self-regulation mechanism: The iron con- 
tent adapts quickly to any changes in the ambient medium (den- 
sity, radiation field) unless the density gets too small out in the 



4 



Peter Woitke: Too little radiation pressure on dust in winds of oxygen-rich AGB stars 




1 2 3 4 5 6 

r/R* 



Fig. 3. Dynamical wind model after 100 years of simulation 
time. Parameter: M* = 1 M Q , = 2500 K, L* = 10000 L Q 
(log g = — 1 .015), Z = 1, piston period P - 600 days and velocity 
amplitude Av = 2km/s. The dotted curves show the hydrostatic 
solution, and the T = 1 level. T is arbitrarily enhanced (see text). 

wind (here at r ^ 5 R+) where the degree of condensation of Fe 
freezes in on a relatively low final level of ~ 17%. It is this factor 
that kills the mass loss (in comparison Mg: ~65%, Al: ~ 100%). 
Noteworthy, silicates in AGB star winds are in fact observed to 
be Fe-poor (Bowey et al. 2002). 

The radiative acceleration just even approaches the gravita- 
tional deceleration around r * 3.5 - 4/?*, which can be con- 
sidered as an analog of the sonic point in stationary winds. 
Determined by long-term averages at the outer boundary, the 
mean mass loss rate, the mean outflow velocity and the mean 
dust-to-gas ratio are (M) ~ 2.3 x 10~ 9 M Q /yr, {v^,) ss 2.6km/s 
and (pd/Pg) * 1.6 X 10~ 3 , respectively (for r = r gas +5 r^st)- 

6. Conclusions 

1. This letter reports on a negative result. Detailed dynami- 
cal models with frequency-dependent Monte Carlo radiative 
transfer and time-dependent formation of dirty dust grains 
cannot explain the observed magnitude of mass loss rates 
from oxygen-rich AGB stars, even in case of extreme stellar 
parameters (i. e. high L*/M* ratios). 

2. The role of solid iron and Fe-rich silicates is crucial for the 
wind driving mechanism. These condensates are the only 
ones that are opaque around 1 fim and, thus, only these con- 



densates can efficiently absorb the stellar light. Since the Fe 
containing condensates are not particularly stable, they form 
at too large distances from the star in order to provide an 
efficient mass loss mechanism. 

3. Previous grey models (Jeong et al. 2003, Ferrarotti & Gail 
2006) have calculated the radiation pressure on dust with 
Rosseland mean opacities which leads to a severe overes- 
timation in the O-rich case, because the mid-IR dust ab- 
sorption resonances are strongly peaked in a spectral region 
where the stellar flux is low. The local flux is not B \{T&)- 
like, because the dust features are optically thick which 
makes the radiation field more isotrop and reduce F^ just 
where Kf u *L is large. 

4. The dust condensation sequence is strongly affected by ra- 
diative transfer effects. Pure, glassy condensates like AI2O3 
have lower dust temperatures than solid Fe or Fe-rich 
silicates. The differences are as large as 1000 K, which 
favours the formation of the glassy condensates and pre- 
vents the formation of Fe-inclusions close to the star. 
The results in this paper are consistent with the obser- 
vational finding of AI2O3 at radial distances as small 
as 1 .5 R* around a Oh (Verhoelst et al. 2006) as well as 
with the observed dust condensation sequence in O-rich 
AGB stars (Blommaert et al. 2006, Lebzelter et al. 2006), 
because AI2O3 can already exist in an extended atmosphere 
without mass loss, whereas the Mg-Fe-silicates form in the 
more distant wind regions which require mass loss. 

5. The mass loss mechanism of oxygen-rich AGB stars and red 
supergiants is still a puzzle. Pulsations alone cannot drive an 
outflow because the radiative cooling of the gas is too ef- 
ficient, even in non-LTE (Woitke 2003, Schirrmacher et al. 
2003). According to the results of this letter, even a com- 
bination of stellar pulsation and radiation pressure on dust 
is insufficient to drive the mass loss. Do we have to re-visit 
Alfven-waves (e. g. Vidotto & Jatenco-Pereira 2006)? 
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